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Abstract 


We study the spherical gravitational collapse of a compact object 
under the approximation that the radial pressure is identically zero, 
and the tangential pressure pg is related to the density p by a linear 
equation of state pg = kp. It turns out that the Einstein equations 
can be reduced to the solution of an integral for the evolution of the 
area radius. We show that for positive k there is a finite region near 
the center which necessarily expands outwards, if collapse begins from 
rest. This region could be surrounded by an inward moving one which 
could collapse to a singularity - any such singularity will necessarily 
be covered by a horizon. For negative k the entire object collapses 
inwards, but any singularities that could arise are not naked. Thus 
the nature of the evolution is very different from that of dust, even 
when k is infinitesimally small. 
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1 Introduction 


The cosmic censorship conjecture remains one of the most important un¬ 
solved problems in classical general relativity. In physical terms the con¬ 
jecture states that singularities arising in generic gravitational collapse are 
hidden behind an event horizon. The lack of a general approach for inves¬ 
tigating censorship has lead to studies of simplihed collapse models, and in 
particular, study of spherical gravitational collapse. The best understood 
spherical collapse is that of dust, where it has been shown that for certain 
initial data the central singularity which forms in collapse is naked, whereas 
for other initial data it is covered |^. A more realistic collapse model would 
invoke pressure. While there are general existence arguments suggesting the 
occurence of naked singularities when pressure is included 0, explicit exam¬ 
ples highlighting the role of initial data are rare. It is thus not known how 
the naked singularity arising in dust collapse is modihed, when pressure is 
introduced. An important exception is the numerical study of self-similar 
collapse of a perfect fluid, where it was shown that naked singularities result 
generically in such a collapse |^. 

The energy-momentum tensor for spherically collapsing matter is diag¬ 
onal in the comoving frame, the components being the energy density, the 
radial pressure, and the tangential pressure. Evidently, the general analyt¬ 
ical solution of the Einstein equations for this system is not known, even 
after an equation of state is specihed. Major simplihcation of the system of 
equations results if the radial pressure is set to zero, and only the tangential 
pressure retained. This is the approximation studied in the present paper. 
Admittedly, it is rather unphysical to consider a collapsing system in which 
there is only tangential pressure. However, our purpose is to explicitly study 
how the introduction of pressure modifies the dust scenario, and the present 
analysis should be considered a small step in that general direction. 

We assume that the tangential pressure pe is related to the energy density 
by a linear equation of state pg = kp where /c is a constant. As a result, 
the Einstein equations can be reduced to a non-linear ordinary differential 
equation for the evolution of the area radius of a fluid element. The initial 
density and velocity prohles of the collapsing object enter as free parameters 
in this equation, which has a one parameter family of solutions, labelled by 
the parameter k. 

The equation can be solved exactly for some specihc values of k, and can 
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be integrated numerically for a general, given, k. However, our aim is to 
check whether any singularities that form are naked, and that can be done 
without resorting to numerical study. The behaviour that we find is quite 
surprising, and quite unlike the evolution of dust. For positive fc, we find 
there is a region near the center which necessarily moves outwards, if started 
from rest, and just cannot collapse. Depending on how far the star extends, 
this outgoing region could be surrounded by a region which collapses inwards, 
and could develop a singularity. However, the singularity will not be naked, 
but covered by a horizon. Beyond this collapsing region, there will in general 
be a second region of expansion. For negative k, the entire object collapses 
inwards, unlike for positive k, but any singularities that arise are covered and 
not naked. Thus the naked singularities in dust collapse are removed when 
a non-zero value of k is switched on, howsoever small. Of course, this may 
not be the case when the radial pressure is included, or when the tangential 
pressure obeys a different equation of state. It does however suggest the 
possibility that when collapse to a naked singularity occurs, the stability of 
the collapse in the presence of small perturbations should be examined. 

In Section 2 we describe the relevant held equations, and in Section 3 
the evolution for positive k (including the case fc = 1/4 which can be solved 
exactly). In Section 4 we describe the evolution for negative k. 

Spherically symmetric static and quasi-static solutions with only tangen¬ 
tial pressure have been studied before, notably by Lemaitre 0, as well as 
by others Those results do not appear to have any direct bearing on or 
connection with the results reported here. 

2 The Field Equations 

In comoving coordinates (t, r, 0, 0) the spherically symmetric line-element is 
given by 


ds^ = — d^dr"^ — B?dVt^ ( 1 ) 

where a and oj are functions of t and r. The area radius R also depends 
on both t and r. In comoving coordinates the energy-momentum tensor for 
a spherically symmetric object takes the diagonal form T/ = {p,Pr,Pe,Pe)- 
The quantities Pr and pg are usually referred to as the radial and tangential 
pressure, respectively. The Einstein held equations for this system are 
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m! = AnpR^R', 


( 2 ) 


rh = —At^Pj-R^ R, 


(3) 


and 


a 


2p; AR' , 

~r-^ t3/ , — \(Ps 

P + Pr R{p + Pr) 


Pr), 


2p _ AR{p + Pe) 
p + Pr R{p + Pr) ’ 


(4) 

(5) 


m = ]^R(l + e-^R^ - e-‘^R'") . (6) 

Here, m{t, r) is a free function arising out of integration of the Einstein 
equations. Its initial value, m(0,r), is interpreted as the mass interior to the 
coordinate r. 

We now make the approximation that the radial pressure is identically 
zero. Eqn. remains unchanged, whereas Eqn. (^) gives the important re¬ 
sult that the mass-function does not change with time, and remains constant 
at its initial value m(0,r) = m(r). This is what simplihes the analysis very 
considerably. Eqn. (|]) remains as above, while Eqns. @ and simplify to 

and 


r[ pj- 


( 8 ) 


Next we assume, for simplicity, that the equation of state is pe = kp, 
where the constant k is chosen to lie in the range — 1 < fc < 1, in accordance 
with the weak energy condition, and the dominant energy condition. Eqn. 
(1^ then has the solution 




( 9 ) 
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where a pure function of time has been set to unity. (This choice merely 
reflects a rescaling of the time coordinate in the metric). Eqn. (|^) has the 
solution 


= V>(r)p2i?4(i+fc). (10) 

Here, 'ip{r) is an aribtrary function of the position coordinate r which we 
retain because we would like to use this freedom later. Using Eqn. (Q) we 
can write this as 




167r2 R'^ ^ 


( 11 ) 


Note that C'(r) is a positive function. Using the solutions (^ and (|^ in 
Eqn. (ID yields the following equation for the evolution of the area radius: 


^2 ^ _ 2m{r)R^^-^ - C{r)R^'^ = 0 . ( 12 ) 


This non-linear differential equation, which can be integrated in principle, 
gives a one-parameter solution of Einstein equations. While the function 
m{r) is determined by the initial density prohle, the function C{r) becomes 
known once the initial velocity prohle is also given. Note that this simplihed 
equation for R results because of the simplifying assumptions that the radial 
pressure is identically zero and that the tangential pressure obeys the linear 
equation of state. If R(t, r) could be solved for from this equation, the 
evolution of the density and the metric components follows from the earlier 
equations. We note that the special case of dust can be recovered from the 
above equation by setting /c = 0. That yields the familiar equation for dust 
collapse. 


(13) 

where we have set C{r) = 1 + f{r). 

We shall now study the properties of Eqn. ([T^ ) and its implications 
for collapse and singularities, treating the cases of positive and negative k 
separately. 
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3 The evolution for positive k 


Evidently, Eqn. (0 cannot be solved analytically for an arbitrary value of 
/c, though numerical integration for a given value of k would be straightfor¬ 
ward. Our interest here, however, is not so much in the exact solution of 
this equation, but rather in those qualitative features which would help us 
determine whether singularities form in collapse, and if so, whether they are 
naked. To this end, we note that for a fixed r Eqn. (^) can be thought of as 
the equation for the one-dimensional motion of a zero-energy particle, along 
coordinate i?, in the potential V(R, r) given by 


V(R, r) = R^^ - 2m(r)R^^-^ - C{r)R^K (14) 

Let us analyse the shape of this potential, as a function of R. By writing 
V{R,r) as 

V{R, r) = R^^-\R -2m- CR^’^+^) = R^’^-^W{R, r) (15) 

we note that V{R, r) will be positive if and only if the function W{R, r) is 
positive. Consider Erst the case of positive k. The extremum of W{R, r) is 
at the point Rq given by 



1 

C'(r)(l + 4A;)' 


(16) 


It is easily checked that this extremum is a maximum. It follows from (p!5|) 
that in order for W{Ro,r) to be positive, the following condition should be 
satished : 


We will assume that collapse begins from rest - it then follows from Eqn. 
0 that C{r) gets hxed as 

Y _ 2m 

C(r) = (18) 

and Eqn. (0 becomes 
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( 19 ) 


/ 2m \ ^ {Akf^ 

[~J (l + 4fc)(i+4^)' 

Note that 0 < 2m/r < 1. It is easily shown that the function {2m/r) 
always satishes this inequality, except when (2m/r) = Ak/{1 + Ak) = fc*. In 
the latter case, the two sides of Eqn. (|T^) are actually equal, and the maxi¬ 
mum value W{RQ,r) is exactly zero. It follows that the maximum W{RQ,r) 
is necessarily positive, except when {2m/r) takes the special value fc*. Fur¬ 
ther, since W{R, r) goes to (—2m) as R goes to zero, and it goes to negative 
inhnity as R goes to inhnity, we conclude that W has two zeroes. These are 
also zeroes of the potential function V{R, r). In the region bounded by these 
two zeroes, V{R,r) will be positive. Thus, as can be seen from Eqn. ([T^) , 
this region will be forbidden, by the requirement of positivity of R^. Outside 
of this forbidden region, both W{R, r) and V{R, r) will be negative - that is 
the allowed region. 

In the region to the left of the lower zero, where V{R, r) is negative, 
the shape of the potential depends on k. For k < 1/4, V{R) monotonically 
decreases as R decreases, ultimately going to —oo as R goes to zero. For 
k = 1/4 the potential again monotonically decreases with decreasing R, 
going to —2m(r) at i? = 0. For k > 1/A the potential will have at least one 
turning point (which is a minimum) and goes to zero at i? = 0. In the region 
to the right of the larger zero, where V{R) is again negative, it will go to 
—cx) as R ^ oo, for all k. While there are no turning points in this region 
for k > 1/4, there may be turning points for k < 1/4, depending on the 
functions C{r) and m{r). The shape of V{R) for the various cases is shown 
schematically in Fig. 1. 

We now show that this shape of the potential function has interesting 
implications for the nature of the motion. A fluid element labelled by r can 
lie either to the right of the larger zero, or to the left of the smaller zero. For 
convenience of analysis, we shall assume that at the time t = 0, which marks 
the start of collapse, we have the scaling R = r. (Recall that such a freedom 
of labelling is available). Since the motion starts from rest, it is evident that 
the potential V{R, r) is zero at the beginning. Hence one of the two zeroes 
of the potential is given exactly hj R = r. It is important to settle whether 
this is the larger or smaller of the two zeroes. For this purpose, we note that 
if the starting point R = r lies to the left of the maximum Rq, it will satisfy 



rP < R^, i.e. 


Ak 2m 

1 + Ak r ’ 


( 20 ) 


as follows after using Eqns. ([T6|) and (^). Similarly, if the starting point lies 
to the right of the maximum, it will satisfy 


Ak 2m 
1 + Ak r 


( 21 ) 



Fig. 1 The schematic shape of the potential function for positive k and for a fixed r. 


Hence for all those points r for which (^) holds R = r is the smaller of the 
two zeroes, and all points for which (|2l|) holds, R = r is the larger of the two 
zeroes. Now consider a neighborhood of the center. The ratio (2m/r) is equal 
to zero at the origin, and hence there will be a neighborhood of the origin 
for which the inequality (^T]) will necessarily hold. All the fluid elements 
in this neighborhood will have their motion begin at the larger of the two 
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zeroes, and such elements will necessarily move away from the origin, towards 
inhnity. They simply cannot collapse. Hence, there will be a neighborhood 
of the origin which will necessarily move outwards, starting from rest. If the 
condition (^T]) is satished all over the star, then all of it will move outwards. 
For fc > 1/4 the shape of the potential implies that everything will escape to 
inhnity. For /c < 1/4 some or all of the huid elements could get trapped at a 
possible minimum the potential could have. 

The ratio 2m(r)/r is zero at the origin, and near r = 0 increases with 
increasing r. Hence it is possible that the above-mentioned neighborhood 
of the center (for which ( pT]) holds) is surrounded by a region in which (p0|) 
holds. If so, this outer region will have the left zero of the potential as 
its starting point. As a result, this region will collapse inwards, and will 
intersect the inner region which is moving outwards. This appears to be 
a shell-crossing singularity. We will assume here that the inner and outer 
regions simply cross each other. In other words, the ingoing region does not 
see the outgoing inner region, effectively. 

This infalling region will necessarily shrink to i? = 0 for /c < 1/4. It could 
do so for fc > 1/4 as well - i.e. it could overshoot the minimum of V{R) for 
suitably chosen velocity and density prohles. i? = 0 possibly corresponds 
to a curvature singularity - this would happen if the energy density dehned 
by Eqn. (2) blows up. In order to actually check whether this happens we 
have to integrate Eqn. (p!2[) to hnd out the behaviour of R' as R goes to 
zero. However, for the present, our intention is not really to check whether 
singularities do form, but to show that if they form, they are not naked. 

We show this by analysing the null geodesic equation for radial geodesics. 
This method has been developed in earlier papers (see, for instance |P), 
and here we briefly recall the equation we need to use. If the coordinate tq 
develops a singularity, the geodesic equation for possible geodesics emerging 
from this singularity can be written as 


X = 


lim 

*ro,R- 


R' 


>0 a\r — 


// + 2m/i?\ 

1 + / )■ 


( 22 ) 


The various new quantities introduced in this equation are dehned as follows. 
X = R/\r — ro|" is dehned as the tangent to a geodesic, and the constant a 
is so chosen that the ratio R '/(r — ro)“~^ is hnite and non-zero in the limit of 
approach to the singularity. The function f{t,r) is dehned using Eqn. (Q) 
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as follows: 


CCrjfl* = !+/((,r). 


(23) 


It can be shown that the singnlarity at r = tq will be naked if this eqnation 
admits at least one positive real root X = Xq. We can easily see that in the 
present case snch a root will not exist. The qnantity R' is non-negative, as a 
resnlt of the weak-energy condition and Eqn. (|]). Note that C'(ro) and m(ro) 
are non-zero qnantities, and that the fnnction / goes to —1 as the singnlarity 
is approached. As a resnlt the expression inside the sqnare-root on the right 
hand side of the above geodesic eqnation goes to oo, thereby preventing the 
occnrence of a positive real root. Hence the singnlarity will not be naked. 

In snmmary, the natnre of the evolntion for positive k is as follows. If the 
evolntion starts from rest, there will decidedly be a region near the center 
which will move ontwards. This is the region in which the ineqnality (PT| ) will 
be satished. In general, this region will be snrronnded by an onter region 
which moves inwards, becanse it satishes the ineqnality (^). This region 
conid become singnlar, bnt snch a singnlarity will not be naked. We note 
two fnrther points. Firstly, the inner and onter regions are separated by an 
infinitesimally thin shell where neither ( ] 20 | ) nor (^) holds, bnt exactly the 
relation 2m/r = 4k/I + 4k is satisfied. For this valne of r, the two zeroes of 
V{R, r) merge and are eqnal to zero, which is also the valne of the maximnm 
of the potential. A particle initially at rest at this valne of r will continne to 
be at rest - this eqnilibrinm is nnstable however. The second point is that the 
entire star conid consist, in principle, of a snccession of ingoing and ongoing 
regions, depending on the ratio 2m/r in these regions. The generality of the 
resnlts derived above will however not be affected, even if this were the case. 

It is instrnctive to compare these results with those for dust collapse, 
where A: = 0. For any value of k, howsoever small, there exists a region 
near the center which will move outwards. In the limit that k goes to zero, 
this region shrinks to zero size. So long as such a region exists, the value 
of m(r) for any infalling shell is non-zero - it is this fact which prevents the 
singularity from being naked. Hence collapse for A: = 0 is qualitatively very 
different from what happens when A: 7 ^ 0. 

We now illustrate our conclusions with the case A: = 1/4, which is exactly 
solvable. The equation of motion ([T^) becomes 


= CR^ + 2m-R 


(24) 
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which has the following parametric solution 


R = 


r — (4m — r) cosh 6 
^ ’ 


t = 


•/c' 


(25) 


We assume that collapse starts from rest. The function C(r) is equal 
to (r — 2m)/r^. The start of the evolution is at t = 0, at which time we 
have R = r. The solution correctly describes evolution for all values of r. 
Points with r > 4m expand outwards, the point r = 4m does not move, 
and points with r < 4m move inwards. There will be a neighborhood of the 
center which satishes r > 4m and this region will move outwards. This is 
consistent with what we expect from the inequalities (1^) and (pT]) above - 
for fc = 1/4 the critical dividing shell between ingoing and outgoing regions 
satishes the relation r = 4m. If the entire star consists only of the region 
satisfying r > 4m then all of it will expand outwards. 

A singularity R = 0 will develop for points with r < 4m at time ts{r) 
given by the singularity curve 


ts{r) = ^^, cosh0,(r) = —(26) 

By calculating R' we can show this is a true curvature singularity, where 
the density and the Ricci scalar blow up. Further, by applying the analysis 
for roots described above, it is shown that this singularity is not naked, but 
covered by a horizon. 


4 The evolution for negative k 


We again analyse the shape of the potential function V{R,r) and of the 

this time for negative k. It 
goes to inhnity. As R 


function fF(R,r) given by Eqns. (lA) and 
can be checked that as R goes to inhnity, W{R, r) 
goes to zero, W{r) goes to —2m(r) for k > —1/4, to —2m — C for k = —1, 
and to —oo for k < —1/4. For k < —1/4, W{R) has no extremum, and for 
k > —1/4 it has an extremum which is a minimum. Thus W{R) will have 
one zero - it will be positive to the right of this zero, and negative to its left. 
This point will be a zero of V(R) as well, and the potential will be positive 
to the right of this zero, thus giving the forbidden region. To the left of this 
zero, V(R) will be negative, hence this is the allowed region. For k < —1/4 
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the potential monotonically decreases to —cx) as R goes to zero, whereas for 
k > —1/4 it could have turning points (including a minimum) before going 
to —CX) at -R = 0. These shapes are schematically shown in Fig. 2. 




Fig. 2 The shape of the potential function for negative k and a fixed r. 


Since the shape of the potential is quite different compared to when k is 
positive, there is no possibility of there being a region which could expand 
outwards, starting from rest. The starting point of evolution is at the zero 
of the potential, which is given hy R = r. From here, the fluid element 
will proceed to a smaller value of R, necessarily hitting i? = 0 for the case 
k < —1/4. For the case /c > —1/4 the collapsing particles will either get 
trapped at a minimum, or hit i? = 0 , depending on the choice of initial 
density and velocity prohles. Thus a curvature singularity could result, for 
r = 0 as well as for r 7 ^ 0. We show now that such a curvature singularity 
cannot be naked. 

For this purpose, we again examine the geodesic equation , using for 
fit,r) the solution given by Eqn. (p3|). Since k is negative, fit,r) goes to 
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infinity as R goes to zero. Hence in the limit of approach to the singularity, 
the expression inside the square-root in Eqn. (|^) is equal to 


2 m(r) 

^ C(r)i?i+4fc' 


(27) 


This quantity is greater than or equal to one, which prevents the roots equa¬ 
tion p^) from having a positive real root. Hence the singularity cannot be 
naked. In contrast, for the case of dust, the quantity inside the square-root 
goes to inhnity for r 7 ^ 0 , but for the singularity at r = 0 , this quantity is 
simply equal to 2m{r)/R (because now /(t,r = 0) = f{r = 0) = 0). De¬ 
pending on how 2m/R behaves in the limit, the r = 0 singularity may or 
may not be naked. 

We now consider the example of fc = —1/4, for which the equation (|T^) 
for the evolution of R can be explicitly solved, and the solution is 


t=-{r- Rf/"^ (2r + R) . (28) 

3 

This solution has been written assuming that collapse begins at f = 0 from 
rest, and the scaling at f = 0 is i? = r. For a given shell r, the singularity 
/? = 0 arises at the coordinate time ts{r) = 4r^/^/3. The quantity R' can be 
calculated, and is equal to (2r — R)/R. Hence from Eqn. (@) we get that the 
energy density evolves as 


m' 

4irfi(2r-ii)' 


(29) 


For r 7 ^ 0, m' is non-zero, implying that at i? = 0 the density blows up, 
and hence there is a curvature singularity. However, the situation is quite 
peculiar for r = 0. At this point, the solution (p 8 | ) is not valid after f = 0, 
which is just the starting epoch of collapse! Hence we cannot talk of the 
density evolution at r = 0 using this solution. 

As for the points r 7 ^ 0, the singularity is not naked. This is seen by 
examining the roots eqn. ( | 2 ^) - the quantity inside the square root is equal 
to {r — R)/{r — 2m), which is greater than one at the singularity R = 0. 
Hence there are no positive roots to this equation. An equivalent but more 
physical way to arrive at this result is to work out the slope dt/dr along a 
light ray and compare it with the slope dts{r)/dr along the singularity curve. 
If the tangent along the light ray has a larger slope than along the singularity 
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curve, it means the singularity lies outside the lightcone, and hence is not 
naked. In the present case, the light ray is given by 

^ )/2 _ ^r-R 

dr y/r — 2m’ 

whereas along the singularity curve we have dts{r)/dr = Thus it 

follows that the slope of the light-ray is larger than that of the singularity 
curve (we recall that 2m(r)/r < 1). 

In this paper, we have not considered the matching to an external Schwarzschild 
solution. If the density and the constant k both vanish outside the boundary 
of the matter, we get the Scwarzschild solution. This would require that k 
actually be a function of r, something that we have not presently considered. 

Our approach for the time being is to consider the behaviour of solutions in 
regions near to the center r = 0, where it should be alright to let k to be 
a constant. In principle, one will have to have to consider matching such 
regions to an exterior matter region where k varies with r, and we leave that 
for a future investigation. 

In conclusion it could be said that in so far as tangential pressure is 
concerned, the naked singularities in dust collapse are very special. They 
disappear when tangential pressure is included. However, these naked singu¬ 
larities will probably reappear when radial pressure is also taken into account. 
Including radial pressure makes the problem analytically untractable, as now 
the mass-function m{t,r) evolves with time. One may have to resort to a 
numerical study to tackle this problem. 

The sign of the pressure appears to play a crucial role in deciding the 
details of the evolution. Similar examples, where the sign of the pressure 
makes an important difference in the nature of the collapse, have been noted 
by Szekeres and Iyer 0 , and also by Cooperstock et al. . Szekeres and Iyer 
found that in order for naked shell-focussing singularities to occur at r > 0, 
it is necessary that the radial or tangential pressure must either be negative 
or equal in magnitude to the density. Cooperstock et al. studied the non¬ 
central shell-focussing singularity for a perfect fluid obeying the weak energy 
condition and showed that with negative pressure, naked singularities can 
form. 

Upon the completion of this work, we received a preprint by Giulio Magli 
[^j who has also examined the dynamics of a spherical matter distribution 
having only tangential pressure. 
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